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Abstract. Boundary value problems for integrable nonlinear evolution PDEs for- 
mulated on the half-line can be analyzed by the unified method introduced by 
one of the authors and used extensively in the literature. The implementation 
of this general method to this particular class of problems yields the solution in 
terms of the unique solution of a matrix Riemann-Hilbert problem formulated in 
the complex fc-plane (the Fourier plane), which has a jump matrix with explicit 
(x, t)-dependence involving four scalar functions of k, called spectral functions. 
Two of these functions depend on the initial data, whereas the other two depend 
on all boundary values. The most difficult step of the new method is the char- 
acterization of the latter two spectral functions in terms of the given initial and 
boundary data, i.e. the elimination of the unknown boundary values. For cer- 
tain boundary conditions, called linearizable, this can be achieved simply using 
algebraic manipulations. Here, we first present an effective characterization of the 
spectral functions in terms of the given initial and boundary data for the general 
case of non-linearizable boundary conditions. This characterization is based on 
the analysis of the so-called global relation and on the introduction of the so-called 
Gelfand-Levitan-Marchenko representations of the eigenfunctions defining the spec- 
tral functions. We then concentrate on the physically significant case of t-periodic 
Dirichlet boundary data. After presenting certain heuristic arguments which sug- 
gest that the Neumann boundary values become periodic as t — >■ oo, we show that 
for the case of the NLS with a sine-wave as Dirichlet data, the asymptotics of the 
Neumann boundary values can be computed explicitly at least up to third order in 
a perturbative expansion and indeed at least up to this order are asymptotically 
periodic. 
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1. Introduction 

This is the second in a series of papers addressing the most difficult problem in the 
analysis of integrable nonlinear evolution PDEs, namely the problem of expressing the 
so-called spectral functions in terms of the given initial and boundary conditions. In 
this problem was analyzed directly, i.e. without employing the so-called Gelfand- 
Levitan-Marchenko (GLM) representations. Here, we revisit the same problem, but 
we employ the GLM representations; it appears that this latter approach might have 
an advantage for studying problems with t-periodic boundary problems. Furthermore, 
for the NLS we also present certain explicit formulas for the large t-behavior of the 
Neumann boundary values corresponding to the simple case of a sine- wave as Dirichlet 
data. 

We refer the interested reader to for an introduction to the unified method of 
[3 [H] and for a discussion of the difference between linearizable versus nonlinearizable 
boundary value problems. Here we only emphasize that the unified method expresses 
the solution q{x, t) of an integrable nonlinear evolution PDE formulated on the half- 
line as an integral in the complex fc-plane. This representation is similar to the rep- 
resentation obtained by the new method for the linearized version of the given PDE, 
but also involves the entries of a certain matrix- valued function M{x,t,k), which is 
the solution of a matrix Riemann-Hilbert (RH) problem. The main advantage of the 
unified method is that this RH problem involves a jump matrix with explicit (x, t)- 
dependence, uniquely defined in terms of four scalar functions called spectral func- 
tions and denoted by {a{k),b{k), A{k), B{k)}. The functions {a{k),b{k)} are defined 
in terms of the initial data qo{x) = q{x,0), whereas the functions {A(k),B{k)} are 
defined in terms of the boundary values. For example, for the nonlinear Schrodinger 
(NLS) equation 

the functions {A{k), B{k)} are defined in terms of go{t) = q{0,t) and gi{t) = qx{0,t). 
However, for a well-posed boundary value problem only one of the boundary values 
(or their combination) is prescribed as a boundary condition; for example, for the 
Dirichlet problem of the NLS, go{t) is given but gi{t) is unknown. The problem of 
determining gi{t) in terms of qo{x) and go{t) is known as the problem of determining 
the Dirichlet to Neumann map. 

For boundary conditions which decay for large t, by utilizing the crucial feature 
of the new method that the dependence of the jump matrix on {x,t) is explicit, 
it is possible to obtain useful asymptotic information about the solution without 
the need to characterize the spectral functions {A(k), B{k)} in terms of the given 
initial and boundary conditions (or equivalently without the need of constructing 
the Dirichlet to Neumann map). However, the complete solution of non-linearizable 
boundary value problems requires the characterization of A{k) and B{k) in terms of 
the given initial and boundary conditions. The effective solution of this problem is 
particularly important for the physically significant case of boundary conditions which 
are periodic in t, since in this case it is not possible to obtain the rigorous form of the 
long time asymptotics, without the full characterization of A{k) and B{k) (in spite 
of this difficulty important results about the large t-asymptotics of such problems are 
presented in |2l |3l i] ) . 

The construction of the Dirichlet to Neumann map can be achieved by analyzing the 
so-called global relation, which is an algebraic equation coupling {a{k), b{k),A{k), B{k)}. 
A breakthrough in the analysis of the global relation was achieved in |5j, where the 
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global relation was solved explicitly for the unknown Neumann boundary values in 
terms of the given initial conditions, the given Dirichlet boundary conditions, and the 
functions {Lj{t, s), Mj{t, s)}'^^^ appearing in the Gelfand-Levitan-Marchenko (GLM) 
representations of the eigenfunctions {^i{t,k),^2{t,k)}. The numerical implemen- 
tation of an improved version of these formulae was presented in [T3], see also |14j . 
Following these important developments, a similar analysis for the sine-Gordon, the 
two versions of the mKdV and the two versions of the KdV, was carried out in [9] 
and [H]. 

1.1. A new formulation of the Dirichlet to Neumann map. Here we first 
present an extension of the above results. In particular, for the NLS we show 
that the global relation actually yields two linear equations coupling {Lj{t, s)}1 and 
{Mj{t, s)}i. The equation for the Neumann boundary values obtained in [5j can be 
obtained by evaluating one of these two equations at s = t. Similar results are valid 
for other integrable nonlinear PDEs like the sG, mKdV, and KdV equations. 

By employing these two new linear equations, it is possible to obtain a simpler 
characterization of the generalized Dirichlet to Neumann map than the one presented 
in [a El [12]. 

1.2. Time-periodic boundary conditions. Boundary conditions which are peri- 
odic in time are of great importance in applications. Examples include the periodic 
motion of a wave maker or the vertical motion of a buoy in the ocean. We will 
consider the following question: Assuming that the Dirichlet boundary conditions 
are periodic in time, do the Neumann boundary values approach a function that is 
periodic in time as t — >• oo? After presenting a heuristic argument why this question 
might be positively answered in the general case, we provide more concrete evidence 
by showing that: (a) For any Dirichlet data go{t) such that gQ{t + tp) — go{t) decays 
faster than as t — t- oo, gi{t) is asymptotically periodic at least to first order in 
the perturbative expansion, (b) For the prototypical example when go{t) is a sine- 
wave, the asymptotics of gi{t) to third order in perturbation theory can be determined 
explicitly and at least up to this order gi is asymptotically periodic. 

1.3. Organization of the paper. In section [2j we derive the new characterization 
of the Dirichlet to Neumann map for the NLS equation. Time-periodic boundary 
conditions are considered in section [Sl 



2. NLS REVISITED 



Let the vector (<I>i, <I>2) satisfy the 'normalized' t-part of the usual Lax pair of the 
NLS evaluated at x = (see [10] for details): 



2k 



go{t) 
Xgoit) 



iX 



\9oit)\' 



-A5i(t) 
\9o{t)\' 



0<t<T, k€C. 

Then, it is shown in [6j that (^1,^2) admits the GLM representations 



^i{t,k) 

<i>2{t,k) 



$1(0, fc) 

$2(0, A;) 



Li{t,s) - -go{t)M2it,s) + kMiit,s) 



ds. 



(2.1) 



(2.2a) 
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^2it,k) = 1 + 



L2it, s) + '-^go{t)Mi{t, s) + kM2it, s) 



(2.2b) 



<t <T, -t <s <t, keC, 



if and only if the functions {Lj(t, s), Mj{t, s)}1 satisfy the following well-posed Gour- 
sat hyperbolic problem: 



Lit - Lis = igi{t)L2 + a{t)Mi + /3(i)M2, 

Mit - Mu = 2go{t)L2 + igi{t)M2, 

L2t + L2s = -iXgi{t)Li - a{t)M2 + XPit)Mi, 

M2t + M2s = 2Xgo{t)Li - iXgi{t)Mi, 0<t<T, -t<s<t, 

Li{t,t)='-gi{t), Mi{t,t) = go{t), L2{t, -t) = M2(t, -t) = 0, 
where the functions a{t) and /3(t) are defined in terms of goit) and gi{t) by 

^{9o{t)giit)-goit)9iit)), /3W = ^ 
The general solution of the equation 

Ft{t,s)-Fs{t,s) = f{t,s), 
is given by (see appendix [B]) 

+ s t + s^ 



= 9 i9o{t)9i{t) - go{t)9i{t)) , m = ^ [m{t) - A|5o(t)|'5o(t)) • 

< t < r, -t<s<t, 



F{t, s) = F ^) + 1^ /(r, t + s- r)dr. 



(2.3a) 
(2.3b) 
(2.3c) 
(2.3d) 

(2.4) 



(2.5) 
(2.6) 

(2.7) 



Thus, equations ( 2.3a| ) and ( 2.3b ) imply that we require the knowledge of the functions 
Li(x,x) and Mi(x,x), x = which are given by the first two equations in (2.4). 
Similarly, the general solution of the equation 

Ft{t,s) + Fs{t,s) = f{t,s), <t <T, -t < s <t, (2.8) 

is given by 

t — s t — 



F{t,s) 



F 



+ / f{T,T + S-t)dT. 



(2.9) 



Thus, equations (2.3c) and (2.3d) imply that we require the knowledge of the functions 
L2{x, —x) and M2{x, —x), x = {t — s)/2, which are given by the last two equations 
in ( [2^41 ). 

In the case of zero initial conditions, s{k) = I, and then the global relation (2.16) 
of section 2 of ^T] yields 

< t < T, G C+, 



<^>i{t,k) = Fi{t,k), ^2{t,k) = F2{t,k), 



(2.10) 



where Fi and F2 are bounded and analytic functions of A; G and of 0{\/k) as 
A: — )• 00 in C^. 

We replace in the first of equations (2.10) ^i by the rhs of (2.2a), we replace s by 
T, we multiply the resulting equation by 

ke^'^^^'-'\ -t<s<t, 

and integrate with respect to dk along dDi which is the oriented boundary of the 
first quadrant of the complex /c-plane: 

ft 



k 



dDi 



-t 



Llit, r) - -go{t)M2it, r) + kMi{t, r) 
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dDi 



ke'^'^''^^~'^Fi{t,k)dk = 0. 



(2.111 



We simplify this equation by using the following facts: (a) -Fi(i, k) is analytic in k in 
C"*" an d of 0(1 /k) as k — )• oo in C^. Thus, Jordan's lemma implies that the rhs of 
(2.11) vanishes, (b) By employing the change of variables 2k'^ = I, we find 



ke'^ik^r-s)^j^ = -6{t - s) 



thus the first two terms on the rhs of (2.11) yield 



vr 
2 



Li{t,s)- -go(t)M2{t,s) 



(c) Integrating by parts the term involving k"^, we find 

rt 



dDi 



e'^'''^^Mi{t,T)dT ] dk 



(2.12) 



1 

2i 



dDi L 



^2ik (t-^)^^(i^ ^) _ e-2»fc'{t+«)Mi(t, -t) 



dk 



The exponential exp[2ifc^(t — s)] is bounded and analytic in Di, thus Jordan's lemma 
implies that the first integral on the rhs of (2.12) vanishes. Similarly, the integral 

ishes. He 

Mlr{t,T) 



with respect to dr from s to t vanishes. Hence, the rhs of (2.12) equals 
'Mi{t, -t 



ICQ 

2 



+ 



-dr 



Co 



By deforming the contour dDi to the real axis we find that 

/vr 



^'^'dl. 



dDi 



Co 



e 4 . 



Hence, equation (2.11) yields 



Li{t,s)-^goit)M2{t,s) + ^ 

^ V ZTT 



Mi(t,-t) 



+ 



Similarly, the second of equations (2.10) yields 

in 



L2(t,s) + -5o(i)Mi(t,s) + 

2 \/2Tr J-t \/s-T 



' Mi,(^,r) 
M2r(t,r) 



0. 



dr = 0, 



(2.13a) 



(2.13b) 



where we have used that M2{t, —t) = 0. 

Evaluating equation ( 2.13b| ) at s = —t and using that L2{t, —t) = 0, it follows that 
Mi{t, -t) = 0. Then, evaluating ( |2.13a[ ) at s = -t, we find that Li{t, -t) = 0. 

In summary, we have derived the following result: 

Proposition 2.1. Let {Lj{t, s), Mj{t, s)}^, < t < T, -t < s < t, be the GLM 

representations associated with the t-part of the Lax pair of the NLS equation (1.1) 
evaluated at x = 0, see equations (2.1)-(2.5). Then, the global relation of the NLS 
yields the following equations: 



Li{t,s)-'-go{t)M2{t,s) + ^ 



Mir{t,T)dT 
\/s - T 



0, 



(2.14a) 
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r / N «^ /N,./ N eT r M2rit,T)dT 

2 V2vr J-t -T 

Mi{t, -t) = Li{t, -t) = 0. 



0, 



(2.14b) 
(2.14c) 



Remark 2.2. Evaluating equation (2.14a) at s = t, we recover the basic formula of 
[5, 1^: 

'■^ Mir{t,T)dT 



gi{t) = go{t)M2it,t) 



e 4 



(2.15) 



Remark 2.3. Equations (2.14) encode the information contained in the global re- 
lation. The advantage of the GLM representations is that the symmetry relations 
associated with the global relation are automatically incorporated into this informa- 
tion. Indeed, letting k — )• —k yields nothing new at the level of GLM representations. 

For a function F{t,s), —t<s<t,we denote by A the Abel transform in the 
second variable: 

F{t,T)^=. 



{AF){t,s)-- 

The inverse A~^ of this transform is given by 



iA-'F)it,s) 



1 d F{t,T)dT 1 



vr ds J_i a/ 



s — T 



IT 



9P{t,T)dT ^ Fit,-t) 



-t 



+ 



S — T 



Vt + s 



Theorem 2.4. Let q{x, t) satisfy the NLS Jj. j[ ) on the half-line and assume that 
qo{x) = 0, X > 0. In this case Mi{t,—t) = M2{t,—t) = 0, thus the inverse Abel 
transform of {Mj}1 is given by 



A-^Mj{t,s) = - f 



'§^it,r)dr ^ 

\/s — T 



j = 1,2, -t<s<t, 0<t<T. 



(2.16) 



The function gi{t) := qx{0,t) is given by 

gi{t) = go{t)M2{t,t) - V2^e-'^'\A-^Mi){t,t) (2.17) 

where the functions {Mj{t, s)}\, \s\ <t < T, satisfy the following system of quadrat- 
ically nonlinear integral equations: 



Mi{t,s)=go 



t + s 



+ 



t + s 
2 



9o{r) (iXgo{r)Mi + V2^e'^/\A~'M2 



+ i(go{r)M2{T,T) - V2TTe-'''/\A~^Mi){T,T)) M- 



dr. 



(2.18a) 



M2it,s) = A 



t -s 
2 



go{r) (igo{r)M2 - ^e'^l\A-^M{) 



(go{T)M2{T,T) - \/2^e*^/^(^-^Mi)(r,r)) Mi 



dr. 



(2.18b) 



Unless otherwise specified, the functions in the integrand of (2.18a) are evaluated at 
(r, t+s — r) and the functions in the integrand of (2.18b) are evaluated at (r, r — t+s). 



Proof. The expression (2.17) for gi{t) follows from (2.15). The first two equations in 
(2.14) can be written as 



Li(t,s) = -go{t)M2{t,s) 



^'-/\A-^Mi){t,s), 
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L2it,s) 



iX 



go{t)Mi{t,s) 



/l-l 



{A-'M2){t,s) 



Using these equations together with (2.17) to ehminate Li, L2, and gi from equations 
(2.3b) and (2.3d), we find the following system for Mi and M2: 

{dt - ds)Mi{t,s) = - 2go{t) (^^go{t)Mi{t,s) + ^e'^/\A-' M2){t, s) 



+ i[go{t)M2{t,t) - V27re-*^/^(^-Wi)(t,t) M2(t,s 



''^/\A-^Mi){t,s) 



iX (go{t)M2{t,t) - V2^e'^/H^''Mi){t,t)) Mi{t,s). 



{dt + ds)M2it, s) = 2Xgoit) ( -go{t)M2it, s) 



The identities (2.7) and (2.9) together with the initial conditions (2.4) yield (2.18). 



□ 



2.1. An effective perturbative scheme. In what follows we show that equations 
(2.17)-(2.18) provide an effective characterization of the Dirichlet to Neumann map, 
namely we show that they yield a well-defined perturbative scheme for computing 
9i{t). 



Substituting into (2.18) the expansions 
Ml = eM^^^ + e^uf^ + • • • , 



Mo 



90 



where e > is a small parameter, we find the following: The terms of 0(e) yield 



The terms of O(e^) yield 

mF'(m) = 4''(^). 



0. 



£\T){A-'Mi\T,T-t + s) 



+ i(^-iMf ^)(r, t)MI''{t, T-t + s 
In general, the terms of 0(e") yield 

M["'\t, s) = ^q"'* ^ ^ ~'~ k^o^^ lower order terms, 
M2^\t, s) = known lower order terms. 



dr. 



Since go{t) is known, it follows that the system (2.18) can be solved perturbatively 
to all orders. Thus, for 'small' boundary conditions, the system yields an effective 
perturbative characterization of Mi and M2, in which each term can be computed 
uniquely in a well-defined recursive manner. 
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3. Time-periodic boundary conditions 

In applications it is common to have boundary conditions which are periodic or 
at least asymptotically periodic in time. For the Dirichlet problem this means that 
go{t + tp) = go{t) for some period tp > 0, or, in the asymptotically periodic case, that 
go{t + tp) — go{t) decays to zero as t — oo. Intuitively we expect the effect of the initial 
data in the boundary region to diminish for large t, so that the leading asymptotics of 
the solution near the boundary should be determined by go{t). This suggests that if 
go{t) is periodic, then the solution q{x,t) will eventually become periodic in time for 
small X as t — 7- oo. In this section, we will investigate the large time asymptotics of the 
Neumann data for the NLS equation, assuming that the Dirichlet data are periodic. In 
particular, we will consider the question whether (asymptotically) periodic Dirichlet 
boundary conditions lead to asymptotically periodic Neumann boundary values, i.e. 
whether the Dirichlet to Neumann map takes (asymptotically) periodic functions to 
asymptotically periodic functions. After presenting a heuristic argument why this 
might be true in the general case, we will provide more concrete evidence in the case 
of vanishing initial data by showing that: (a) For any Dirichlet data go{t) such that 
go{t + tp) — go{t) decays faster than l/\/t as t ^ oo, gi{t) is asymptotically periodic 
at least to first order in perturbation theory, (b) For the special case when go{t) is 
a sine- wave, the asymptotics of gi{t) to third order in perturbation theory can be 
determined explicitly and at least up to third order gi is asymptotically periodic. 

3.1. Heuristic argument. In order to establish the asymptotic periodicity of gi{t), 
we must analyze the analogous property of the solution {Mi{t, s), M2{t, s)} of equa- 
tions ( 2.18| ). These equations involve the following three types of integrals: 

{IiM){t,t + s):= [ g{T)M{T,t + s-T)dT, (3.1) 

Jt±s 



{hM){t,t-s):= I g{T)M{T,s-t + T)dT, (3.2) 

(/3M)(t, s) := r ^^0^dT, 0<t<T, -t<s<t. (3.3) 
J-t Vs-T 



We assume that g{t) is 27r periodic and analyze the properties of the integrals (3.1) 



(3.3), under the transformations 

t^t + 27r, s + 2tt. (3.4) 



We start with the integral in (3.1): 

r-t+2TT 



{IiM){t + 2Tr,t + s + Att) = g{T)M{T,t + s + Att - T)dT. (3.5) 



Replacing r with t + In in the rhs of the above equation we find 

(/iM)(f + 27r,t + s + 47r) = {IiM){t,t + s), (3.6) 

where 

M{t,s) = M{t + 2TT,s + 2TT). (3.7) 



Similarly, replacing r with r -|- 27r in the integral resulting from (3.2) under the 



transformations (3.4), we find 



{l2M){t + 2TT,t- S) = [ g{T)M{T + 2TT,S-t + T + 27r)dT 
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[ g{T)M{T + 2-K,s-t + T + 27r)dr 



+ 



t-s 
2 



g{T)M{T + 2'K,S-t + T + 2TT)dT. 



-2lT 



Replacing in the second integral in the rhs of the above equation r with (t — s)/2 + r, 
we find 



{I2M) (t + 2TT,t - S) = il2M){t,t - S) 



(3.8) 



+ 



2tt 



+ t]M 



+ r + 27r, 



t 



+ r + 27r )dr. 



We note that for large t the function M appearing in the integral of the second term 
of (3.8) can be approximated by M(t, —t) and also recall the equations 

Li{t, -t) = Mi{t, -t) = M2{t, -t) = 0. 

Finally, we consider the integral in (|3.3l). Replacing r with t -\- 2tt in the integral 



resulting from (3.3) under the transformations (3.4), we find 



(/3Af)(t + 27r,s + 27r) = J 



Mr{t + 2TT,T + 2TT) 



-An Vs-T 

Mr{t + 27r,r + 2tt) 



dr 



dr + 



Mr{t + 27r,r + 2tt) 

\/s — T 



dr. 



-t VS -T J-t-iTT 

Replacing in the second integral in the rhs of the above equation r with r — t, we find 

''^ Mr{t + 27r,-t + T + 2tt) 



{hM){t + 27r, s + 2tt) = {hM){t, s) + 



47r 



^/s + t■ 



-dr. 



(3.9) 



We note that if is bounded in t, the integral in the rhs of (3.9) vanishes as t — )• 00. 

The above properties suggest that it might be possible to prove directly from 
equations (2.17)-(2.18) that the periodicity oi §^{1) implies the asymptotic periodicity 
of gi{t). However, the rigorous proof of this result remains open. 

3.2. Time-periodic boundary conditions in the linear limit. 

Theore m 3.1 . Let q{x,t) = eqi{x,t) + 0(e^) be the perturbative solution of the NLS 
equation (1.1) on the half-line with g(x, 0) = 0, > 0, and q{0,t) = egQi{t) + O(e^), 
where goi{t) is a sufficiently smooth function compatible with the zero initial data, i.e. 
goi and its derivatives vanish att = 0. Let gi{t) = qx{Q,f) = egii{t) + 0{e^), denote 
the Neumann boundary values. If 

goi{t + tp)-goi{t) = o( ^ 



t2- 



t — )• 00, 



for some tp > and 1^ > 0, then 



gii{t + tp) - gn{t) = < 



Oit-"), 0<z^<l/2, 
0(^-1/2 Int), u = l/2, 
[0(ri/2), ^>i/2, 



00. 



(3.10) 



(3.11) 



Proof. The Neumann boundary values gii{t) are given by an Abel transform of goi{t), 
see ^11] : 



9ii{t) 



e 4 



/i(t), where h{t) :- 



m{s) 



ds. 



(3.12) 
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We write 

h{t + tp) 



t+t 



p goi{s)ds 



goi{s)ds 
-\- tp — < 



+ 



t+t 



" goi{s)ds 
-y/t + — . 



Changing variables s — )• s + tp in the second integral, we find 



h{t + tp) 



goi{s)ds 
\/t + tp — s 
goi{s)ds 

\/t + tp — 5 



+ 



501 (s + tp)ds 



Vt- s 







Since the first integral on the rhs is of 0{t ^^'^) as t — )• oo, we infer that 

501 (s + tp) -5oi(s) 



h{t + tp) - h{t) 







-ds + 0{t' 



-1/2N 



t — )• oo. 



Suppose now that 501 (t) satisfies (3.10). Then, there exist constants K > {) and 
M > 0, such that 



* goi (g + tp) -goi(^) 
Vt-s 



ds 



< 



K 



\gm{s + tp) -goi(s)| 
Vt- s 



/ ^ ds 

Jk 82+" Jt - S 



10 V--*' JKsi^^^/T-S 

for t > K. The first integral on the rhs is clearly of 0(t~^/^) as t — t- 00, whereas the 
second integral satisfies 

ds 



K s^+^y/ir^ 



V2 2' 2' t J t^ 

'0(t-^), 0<i/<l/2, 
0(t-V2int), z^ = l/2, t- 



00. 



This proves (3.11 ). 



□ 



3.3. A sine- wave as boundary data. The prototypical example of a periodic 
boundary profile which vanishes at t = is the sine-wave 5o(t) = sint. For this 
example we can take the results of the previous subsection one step further: We 
will compute the asymptotics of gi{t) to third order in perturbation theory and, as 
expected, find that the result is asymptotically periodic. 



Theorem 3.2. Let 



q{x, t) = €qi{x, t) + €^q2ix, t) H 



0, 



be the perturbative solution of NLS equation (1.1) on the half-line with the initial 

(3.13) 



conditions q{x, 0) = 0, x > 0, and the Dirichlet boundary conditions 

q{0,t) = esint + 0(e''). 
Let 

5i(t) = Q^(0,t) = e5ii(t) + eV2(t) + eV3(t) + 0(e^), 
denote the corresponding Neumann boundary values. Then 



9ii{t) 



e 4 



(cost + sint) + 0(t 



-3/2x 



5l2(t) = 0, 

'^5'13(0 = ci cost + C2 sint + C3(cos3t — sin3t) + 0{t^^^^^) 



t — )■ 00, 

t > 0, 
t — )• 00, 



(3.14a) 

(3.14b) 
(3.14c) 
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where the constant C3 is given by 



C3 



:l-^/3)(^/3 + - 
16 



whereas expressions for the constants ci and C2 are given in equation (3.43) below. 

Equations (3.14) imply that the Neumann boundary values are asymptotically pe- 
riodic at least to third order in the perturbative expansion. More precisely, 



5ii(t + 27r)-gn(t) = 0(t-3/2) 
5i2(t + 27r)-gi2(t) = 0, 



t — )• 00, 
t > 



and 



00. 



(3.15a) 
(3.15b) 

(3.15c) 



gi3(t + 2^)-5i3(i) = 0(i-'/''), t 
Remark 3.3. If higher order terms are included in ( |3.13 ), i.e. 

q{0,t) = esint + e^go2it) + e^gosit) + O(e^), 

where (702 (i) and gosit) are sufficiently smooth functions compatible with the zero 
initial data, then the only effect of the functions go2{t) and go3{t) is to add the terms 



* 502(5) 



ds and 



* 503 (s) 



ds 



to the expressions for 512 (i) and 513 (t), respectively. Employing the argument used 
to prove theorem 3.1, we deduce that gi{t) is still asymptotically periodic to third 
order in the perturbative expansion provided that 502 (0 and 503 (i) are asymptotically 
periodic. 



Proof of theorem 3.2. Since A only affects the sign of 513, we may assume that A = 1. 
When goi{t) = sint, the function h{t) defined in (3.12), i.e. 



hit) 



cos s 



Vt 



zds, 



satisfies 



h{t) = 2ViiF2 (i;^,^; 



(3.16) 



^(sint + cost)-^ 



+ 



15 



945 



8t7/2 32tii/2 



+ 0{f 



-15/2N 



00. 



Since gn = — e~T/i/y^, this proves (3.14a). Since goi{t) has no term of O(e^), gi2{t) 
vanishes identically, which proves (3.14b). Clearly, equations (3.15) are a consequence 
of (3.14). Therefore, it only remains to prove (3.14c), and this relies on the follow- 
ing expression for 5113 (t), which is obtained by setting goi{t) = sint in the general 
expression for 513 (t) derived in [TT] : 

cost" 



A5i3(i) 



2c 
iri 



* sin^f 
Vt^ 



dt' 



2c 
vri 



■ 2 ./ 

sm t 



dt"dt' 



— \ — sint [ h{t')sint'dt' 
V vr Jo 



c c ct c 

— Issc{t) + -^hhhit) + -^hshit) + -^Ihhs{t) 

TTi vr^ vr^ vr^ 



(3.17) 



=: Ti + r2 + Ts + r4 + r5 + Te + Tj, 
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where c = 4 and the functions Igsc, Ishh-, Ihsh: and Ihhs are defined by 



Ishh{t) 
Ihsh{t) 



t' 



dt' sin t' I dt" sin t" 



cos r 



{t-f + t" - 

fdt'sint' f dt"h{t") f dt'" ^^^'"^ 
Jo Jo Jo 



dt'h{t' 



t' 



dt sm t 



dt" 



Ihhs 



ft ft' ft" 

(t) = / dt'h{t') / dt"h{t") / dt" 

Jo Jo Jo 



{t-f + 1" - t"'fl^ ' 

h{t"') 

{t-f + t" - ' 

sin(t"0 
{t-f + t" - t'")3/2 ' 



(3.18a) 
(3.18b) 
(3.18c) 
(3.18d) 



We win treat the different expressions in (3.18) in a series of lemmas. In this respect, 
we introduce the foUowing notations: 



S{z) and C{z) are the Fresnel integrals defined by 



S{z) 



TTS 

sm ds, 

2 



C{z) 



cos —^ds. 



Given a function f{t), Af{t) denotes the Abel transform, 

Am :- ^W'*' 



/o \/t - s 

Given a function f{t), I f{t) denotes the integral 



If{t) ■■= / f{s)ds. 
Jo 

The functions Hs{t) and Hc{t) are defined by 

Hs{t) := I h{s) sin sds, Hc{t) := / h{s) cos sds. 
Jo Jo 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



Note that the inversion of the Abel transform yields Ah{t) = irsmt. The function 
h{t) can be written as 

h{t) = V2^(^c(^^Vi^ cos(t) + sin(t)^ . (3.23) 

Lemma 3.4. The term Ti, i.e. the first term of the rhs of (3.11), satisfies 

Ti 



3V27r 



-^c{J^Vt^ sin(t) + V2,c(^^j^Vt^ sin(3t) 
+ 9s[J^Vt^ cos(t) - V^s(^^j^Vt^ cos(3t)|. 



Proof. 



* sin^ f 1 /■* 3 sin f - sin 3t' , 1 /"* 3 sin(t - s) - sin(3t - 3s) 



, --dt' ^ 

Vt^ 4 



Vt^ 



-dt' 



-ds. 



The result follows by expanding the sine functions and expressing the result in terms 
of C(z) and 5(z). □ 
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Lemma 3.5. The terms T2 and T3, i.e. the second and third terms of the rhs of 

t sin 2t ' 



{3.17 ), satisfy 



To 



2c 
vri 



, / , . f t sin 2t 

,h- A[ cos(t) 

2 4 ; V V2 4 



and T3 = —y j^Hs{t) sinf. 

Proof. This follows using integration by parts and the identity 

2 , ^1 

sin sds = 7 sin 2t. 

2 4 



□ 



In order to find the asymptotics of T^-Tj, our strategy is to reduce the expressions 
for Issc, Ishh, Ihsh, hihs to expressions involving only the Fresnel integrals C, 5*, the 
functions Hc^ Hg, h, and the operators A and /. 



Lemma 3.6. The function Issc{t) defined in (3.18a) can he expressed in the following 
form: 

+ \/27rC ^ (3 sin t + sin(3t ) - 2t cos t) 

+ 2 sin t (^y2TTS + sin 2t) + 2\^t cos | . 

Proof. The relevant integrals can be computed using the type of argument employed 



in the proof of lemma 3.4 and integration by parts 



□ 



Lemma 3.7. The function Ishh{t) defined in (3.18b) can be expressed in the following 
form: 

Ishh = 2A{Hsh) + 2cos{t)I[hA{HcCos)] + 2 sm{t)I[hA{H s cos)] 
+ 2 cos{t)I[hA{Hs sin)] - 2 sm{t)I[hA{Hc sin)] 

- 2A[cos{t) I {hHc cos)] - 2A[sm{t) I {hHs cos)] 

- 2A[cos{t)I{hHs sin)] + 2A[sm{t) I {HHc sin)] 

- 2HcA{Hc cos) - 2HcA{Hs sin) + 2A{hI cos) 

- 2HsA{Hs cos) + 2HsA{Hc sin) + 2A{H'1 cos). 



(3.24) 



Proof. Exchanging the order of integration in the definition (3.18b) of Ishh and then 
computing the integral with respect to dt' , we find 



shh 



* /■*" /■* sini' 

dt"h{t") / dt'"h{t"') / dt'- — — 

Jo Jt" {t-f + t" -t"y/^ 



dt"h{t" 



dt"'h[t"')[4^-:^^ 

^X^t" -t'" 



(3.25) 



+ 2V27rcos(t + t" -t" 



C\ yl-Vtr~^' 



c{ \I-Vt^' 

IT 
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+ 2V27r sm{t + t" - t'' 



Integration by parts implies that the terms involving C ( y ^ \/t" — t'"^ and S{^\J y/W 
can be written as 

2 f dt"h{t") f dt"'( f /i(r) cos(t + t" - r)dr^ " 
Jo Jo \Jo 



and 



t ft" / 

II I. 1 4.11 \ I Jj./// / 



yjt" - t'" 

2j^dfh{t")j^ df'^J^ h{r) sin(t + t" - r)dr^ "'""J^r^^ , 

respectively. Using the identity 

cos(t + t" - r) cos(t" - t'") + sm{t + t" - r) sm{t" - t'") = cos(r - t - t'"), 
it follows that the total contribution of these two terms is given by 

2 f dt"h{t") f dt'" ( f h{r) cos(r - t - t"')dr^ ^ 
Jo Jo \Jo 



Vt" - 1" 



Similarly, integrating by parts twice, we infer that the terms on the rhs of (3.25) 
involving C[^J^\/t — t'") and S ( \/^V t — t'") can be written as 

cos(t - t") 



t / ft" 



and 



2 /* dt"h{t") f dt'" ( f h{r) cos(i + t" - r)dr\ "^" ^ ~ ^'" ^ 
Jo Jo \Jo J Vt — t'" 



2 dt" h{T) h{r) sin(t + r - r)dr^ dr^ ''""^^ ^"^ 



t ft" / ft" 

II 1 I ji\ I JJ.III I 



Vt^ 

-2 J dt"h{t") j dt"'^j h{r)s\n{t + t" -r)dry^^^Z^, 

respectively. Using the identities 

cos(t + T — r) cos(t — t") + sin(t + r — r) sin(t — t") = cos(r — t" — r), 
cos(t + t" — r) cos(t — t'") + sin(t + t" — r) sin(t — t"') = cos(r — t" — t'' 

we infer that the total contribution of these two terms is given by 

- ' 11 ill" "'^^ {[ - - ^'*-) Tfb 



2 fdt"h{t") f dt"'( f 
Jo Jo \Jo 



h{r) cos(r — t" — t"')dr 



In summary, 



' dt"h{t") f dt'"h{t'")^ ^''""^ ^ 
Jo 



^t" - 1'" Vt^' 



(3.26) 



+ 2 / dt"h{t") [ dt'" I [ h{r) cos(r - t - t"')dr 
Jo Jo \Jo 

-2j^dt"(^J^ Hr)(^J^ /i(r)cos(r-t"-r)(ir^(iT 



1 



Vt" - 1'" 
1 

Vt^ 
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,11 J /,//\ / 7,/// I 



1 



2 1 dt"h{t") / dt"'[ / h{r)cos{r-t" -t"')dr 

Jo Jo \Jo ) yt — t" 

The term involving J'^f™l„i equals 

2sm.{t)I{hAh) = 27rsin(i)F^. 
Integrating by parts with respect to dt" , we can write the term involving — as 

-2Hs{t){Ah){t) + 2{A{Hsh)){t) = -2^Hs{t) sin(t) + 2{A{H,h)){t). 
Moreover, using that 

cos(r — t —t ) = cost cos(r — t j+smt sm(r — t j 



and integrating by parts with respect to dt" , we can write the last line of (3.26) in 
the following form: 

rt / rt'" 



h{r) cos(r — t"')dr 



1 

Vt^' 



-2H,{t) l^dt"'(^l^ 

+ 2 / dt"HJt") / h(r) cos(r - t")dr , 

Jo Jo Vt^ 

-2Hs{t) dt'"(^J^ 



4-' 

h{r) sin(r — t"')dr 



+ 2 [ dt"Hs{t") [ h{r)sm{r -t")dr ^ 
Jo Jo vt — t" 

- 2Hc{t)A{HcCos){t) - 2Hcit)A{Hs sin)(t) + 2A{H^ cos)(t) + 2A{HcHs sin)(t) 

- 2Hs{t)A{HsCos){t) + 2Hs{t)A{HcSm){t) + 2 A{H] cos) (t) - 2A{HsHcSm){t). 



Expanding the cosine functions in the second and third lines of (3.26) in a similar 



way, we find (3.24). 



□ 



Lemma 3.8. The function Ihsh{t) defined in (3.18c) can be expressed in the following 
form: 

= AA{hHs) - 2ttHs sint - 2,A{HI cos) - 6A{HcHs sin) + 3^(ij2 cos) 
+ 2A[cos{t) I {hHc cos)] + 2A[sm{t)I{hHs cos)] 

- 2A[cos{t)I{hHs sin)] + 2A{sm{t)I{hHc sin)] (3.27) 
+ 2HcA{Hc cos) + 2HcA{Hs sin) + 2HsA{Hc sin) - 2HsA{Hs cos). 



Proof. Exchanging the order of integration in the definition (3.18c) of Ihsh and then 
computing the integral with respect to dt" , we find 

hsh= f dt'h{t') f dt'"h{t'") f dt"- 



JO 
t ft' 
dt'h{t') / dt'"h{t'") 



sin t" 



t'" {t-t' + t" - t'")3/2 

2sint'" 2sint' 



yjt - t' - t" 



+ 2 V27r / dt'h{t') f dt"'h{t"') cos(t -t' - t" 
Jo Jo 



X ( C( \I-Vt^t 

TT 



"' 1 - c( \i-Vt^' 

vr 
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+ 2^ fdt'h{t') f dt"'h{t"')sm{t-t' -t'") 
Jo Jo 



Exchanging the order of integration, integrating by parts with respect to dt' , and then 
changing back the order of integration, we find that the terms involving C(^^J^\/t — t'^ 

and <S'(-i/ ^^/t — t'^ can be written as 



-2 dt' dt"'h{t"') U^' h{r) cos(i - r - t"')dr 



cos(t - t') 



and 

sin(t-t') 



-2 / dt' [ dt"'h{t"') { [ h{r) sin(t 
Jo Jo \ Jt'" 



r - t"')dr 



Vt^ ' 



respectively. Using the identity 

cos(t - r - t'") cos(t - t') + sm{t - r - t'") sm{t - t') = cos(r -t! + t'"), 
it follows that the total contribution of these two terms is given by 

-2 J -|L= j dt"'h{t"') drh{r)cos{r-t' + t"'). 

Similarly, integrating by parts twice, we infer that the terms involving C(y^\/i ~ t" 
and S{^J^\/t — t'") can be written as 

2 ^* dt' ^* h{T) i^j^^ h{r) cos(i - r - r)dr j dr^ ~ 



^/t^ 

+ 2 /* di'" ( f h{r) cos(t - - r)dr^ """'^^ " 

Jo \Jo 



^/t-t' 



III 



and 



sin(i - t') 



Vt^ 

^* ' ^ sin(t - t'") 



ft ft / ft 

+ 2 dt'h{t') / dt"'( / h{r)sui{t-t' -r)di 
Jo Jo \Jq 



Vt^^ ' 

respectively. Using the identities 

cos(t — r — r) cos(t — t') + sin(t — r — r) sin(t — t') = cos(r — t' + r), 
cos(t -t' -r) cos(t - t'") + sin(t - t' - r) sin(t - t'") = cos(t' - i'" + r) 

we infer that the total contribution of these two terms is 

+ 2 ^* dt'" (^J^ h{r) cos(t' - t'" + r)cir^ ^ ^^^^ . 
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In summary, 



Jo JO 



(3.28) 



2 /* dt' f dt"'h{t"') ( f 
Jo Jo \Jt"' 



h{r) cos(r - t' + t"')dr 



1 



+ 2 / dt'h{t') / dt"'( / /i(r) cos(t' - t'" + r)dr 
Jo Jo \Jo 



1 

1 



Integration by parts with respect to dt' implies that the term involving , equals 



-2Hs{t)Ah{t) + 2A{hHs){t) = -2i:Hs{t) sint + 2A{hH,){t). 



We next consider the second line of (3.28): the integrand in the factor 

rt' ft' 

dt"'h{t"') / drh{r)cos{r-t' + t"') 
10 Jf" 

is symmetric in r and t'" , therefore we can write this factor as 

t' rt' 



1 

'0 Jo 



dt'" / drh{t"')h{r) cos(r - t' + t'"). 



Expanding cos(r — t +t ), we see that the second line can be written as 

-A{Hl cos) - 2A{HcHs sin) + A{H^ cos). 
The third line can be handled by expanding cos(r — t' + r). Finally, using the identity 
cos(t' — t'" + r) = cost' cos(t"' — r) + sint' sin(t"' — r) 



and integrating by parts with respect to dt' , we can write the last line of (3.28) in 
the following form: 



2i/e(t) r^£^ r drh{r) cosit'" -r) -2 f dt'Hcit') f drhirf''''^*' 
Jo vt — t'" Jo Jo Jo 



+ 2Hs{t) 



, drhir) sm{t"' - r) - 2 I dt'Hs{t') I drh{r f''fzl~^ 

yjt — t' Jo Jo Jo yjt — t 



t ^^m rt 



t i-t' 



= 2HcA{Hc cos) + 2HcA{Hs sin) - 2A{H^ cos) - 2A{HcHs sin) 
+ 2HsA{Hc sin) - 2H,A{H, cos) - 2A{HcHs sin) + 2A{H'^, cos). 



Summing up the various contributions, we find (3.27) 



□ 



Lemma 3.9. The function Ihhs{t) defined in (3.18d) can be expressed in the following 
form: 

Ihhs = 2A{hHs) + 2A{hI cos) - 2A[cos{t)I{hHc cos)] - 2A[sm{t)I{hHs cos)] 
+ 2A{hI cos) - 2A[cos{t)I{hHs sin)] + 2A[sin(t)I(/ii?c sin)] 



cos{t)I{hHc cos) + sm{t)I{hHs cos) 
— cos{t)I(hHs sin) — s'm{t) I (hHc sin) 
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+ 2V2^Ts(^^J^^/i^ - cos{t)I{hHs cos) - sm{t) I {hHc cos) 



+ 2 / dt'h{t') / dt 



+ cos{t) I {hHc sin) — sm{t)I{hHs sin) 
,,cost"Hc{t") +smt"Hs{t") 



Vt-t' + 1" 



(3.29) 



Proof. Computing the integral with respect to dt'" in the definition (3.18d) of Ihhsi 
we find 



hhs = / dt'h{t') / dt"h{t") 



t' 



2sint" 



+ 2V2^ [ dt'h{t') I dt"h{t") cos(t -t' + t" 
Jo Jo 



c \ -Vt^] -c \ -Vt-t' + t" 



+ 2V2^ [ dt'h{t') [ dt"h{t") sm{t -t' + t" 
Jo Jo 



s \ -Vt^] -s \ -Vt-t' + t" 



Exchanging the order of integration, integrating by parts with respect to dt' , and then 
changing back the order of integration, we find that the terms involving C(y^J^\/t — t') 

and S{\ -\/t — t') , can be written as 



2 [ dt' [ dt"h{t") ( [ h{r) cos(t - r + t")dr 
Jo Jo \Jt" 



and 



t rt' 



t' 



2 dt' dt"h{t") ( / h{r) sin(t - r + t")dr 
Jo Jo 



cos(t - 1') 

Vt^ 

sm{t - t') 
Vt~F ' 



respectively. Using the identity 

cos(t - r + t") cos(t - t') + sin(t - r + t") sin(t - t') = cos(r - t' - t"), 
we see that the total contribution of these two terms is 



2 I dt' j dt"h{t")( I h{r)cos{r-t' -t")d^^^. 
Jo Jo \Jt" J yt-t' 

Similarly, integrating by parts with respect to dt" , the terms involving C[\/ ^Vt — t' + 

and S(y^J~^^/t — t' + t") can be written as 

- 2V2^c(^\j^Vt^ dt'h{t') h{r) cos(t - t' + r)dr 



+ 2 / dt'h{t') I dt"[ / h{r)cos{t-t' + r)dr 



and 



- 2V2^s(^^Vt^ dt'h{t') h{r) sm{t - t' + 



cos(t -t' + 1" 
Vt-t' + 1'^ 

r)dr 



THE UNIFIED METHOD II 



19 



+ 2 / dt'h{t') / dt"( / h{r) sm{t - t' + r)dr 



sm{t -t' + 1") 
~7t-t' + t" 
Using the identity 

cos(i — t' + r) cos(t — t' + t") + sin(f — t' + r) sin(t — t' + t") = cos(r — t"), 
we infer that the total contribution of these two terms is 

- 2V2^c(^yj^Vt^ dt'h{t') h{r) cos(t - t' + r)dr 

- 2V2^s(^^J^Vi^ dt'h{t') h{r) sm{t - t' + r)dr 



In summary, 

hhs = I dt'h{t') I dt"h{t") 
H rt' 



2smt" 

t' 

h{r) cos(r — t' — t")dr 



+ 2 f dt' f dt"h{t") ( I 
Jo Jo \Jt' 

- 2V2^c(^^Vi^ dt'h{t' 

- 2V2^S ^) ^ dt'h{t') h{r) sin(t - t' + r)dr 



) / h{r) cos(t — t' + r)dr 



+ 2 I dt'h{t') / dt"[ / h{r)cos{r -t")dr 



1 



^t-t' + 1" 

The first term on the rhs equals 2A{hHs). Expanding the trigonometric functions in 
the remaining terms, we arrive at (3.29). □ 

We will find the asymptotics of 513 by determining the asymptotics of the various 
terms appearing in the above lemmas. We begin with the functions He and Hg, which 
can be expressed in terms of C and S. 



Lemma 3.10. The functions He and Hg defined in (3.22) satisfy the following esti- 
mates: 

Hc{t) = ^|\/2^C^y|^\/t^(2t + sin2t) - V2^s(^^^Vt^ cos2t - 2\/tsin(t) 

1 pK 1 /vr , . „ „ , sin(t) 3cos(t) 15sin(t) ^, q/n, 

= 2 V 2' + 4 V 2 - <=»20 - ^ + ^ + + 0(r»/^), 

Hg{t) = ^|2\/^cos(t) - \/2^^C^y^\/t^ cos(2t) + s[\j^Vt^ (sin2t - 2t) 

1 pK 1 /vr, , cos(f) Ssinft) 15cos(t) ^, n/n. 



Proof. This is a direct consequence of using the expression (3.23) for h{t). 



□ 



Since He and H3 are of 0{t) as t — t- 00, by scrutinizing the above lemmas it becomes 
clear that in order to determine the asymptotics of 513 to order 0(t~^/^^), we need 
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to find the asymptotics of the following terms to 0(t~^/^^): 

A{Hsh), A[cos{t)I{hHcCos)], A[sm{t)I{hHsCos)], A[cos{t) I {hHs sin)], 
A[sm{t)I{hHcSm)], A{H^ cos), A{H^ cos), A{HcHsSm), 
I[hA{HcCos)], I[hA{Hs cos)], I[hA{Hs sin)], I[hA{HcSin)], 
I{hHcCOs), I{hHsCOs), /(/liJ^sin), I{hHcSin), (3.30) 

as well as of the following terms to 0{t~^'^/^^): 

A{Hc cos), A{HsSin), A{HcSin), A{Hs cos). (3.31) 

The following lemma will be the main tool for finding the asymptotics of the relevant 
Abel integrals. 

Lemma 3.11. Let N > 0. Let f{t), t > 0, be a smooth function which is bounded as 
t 0, and which satisfies 

f{t) = fa{t) + 0{t-^), i^OO, 

where 

f n.\ ^dff , V^^ +1 • I /i + En=i(cinCosnt + si„sinnt) 
= 1 /?■ + Z^(cj" cosnt + Sjn smnt) | + , 



j=0 \ n=l 



and fj, Cjn, Sjn, cin, Sin, j = 0, . . . , 4, u = 1, . . . , N , and fi are constants. Then, 
the asymptotic behavior of Af(t) as t oo is given by 

Am = zj:^, ^ + M 



2 

4 N 



/— 1 iv ^ 
I E E - 2)^- - ^j-t + 32nH^) sinnt 

j=0 n=l ^ 

+ (-3(j - 2)j + 8jnt + S2nH'^) cos ntj 

— 1 i\ 

I E E 7%*^~'(0' + 4^*) ^^^^ + - 4"*) + o(t-'/''). 

j=0 n=l 

Proof. Since /(t) = O(t^) as t ^ oo, there exist constants K > and M > such 
that 



4 AT 



-4^ds < r ^^ds + M r" -^ds, > (3.32) 

7o yt-s Jo Vt-s Jk Vt-s 

The first integral on the rhs is of 0(t~^/^) as t ^ oo, while the second can be 
computed explicitly and is of 0{t~^^^) as t — t- oo. 

On the other hand, by assumption, there exists L > such that \f{t) — /a(i)| < j 
for t large enough. Thus, as t ^ oo, 

Jt^/a y/t-S Jt^/»Sy/t-S 



"K \ f/^\\ i-t^/^ ^2 
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We conclude that to 0{t ^^^) the asymptotics of Af coincide with the asymptotics 
of J^i/g ^j==ds. In view of the estimates 



smns 



^1/8 y/SyJt — S 

cos ns 



tl/S y/syjt — S 



< 



< 



smns 



ds = o{t-^^^^), 



cos ns 



ds = 0{r 



^l/16^ 



the contribution to the asymptotics of the terms involving cin and si„ are of 0(t~^/^^). 
The contribution of the term involving /i is /ivr + 0(t~^/^^). The contributions of 
the remaining terms can be obtained by expressing the integral^ 



/ { fi ( Cjn COS nt + Snn sin nt) 
Jo \' J 



i = 0,...,4, 

in terms of hypergeometric functions and computing the asymptotics. 



□ 



Although lemma 3.11 is sufficient for most purposes, we need a more sophisticated 
argument to deal with the higher order asymptotics of the terms in (3.31). 



-(sint + cost) 



vr sin t it cos 3t 



Lemma 3.12. Let A,Hc,Hs be defined in {3.20^ and {3.2^ . Then, 

TTf 

A{H,cosm = - 
■jrf 

A{HsSm){t) 
A{H,sm){t) : 
A{HsCos){t) 



vrsint 7rcos3t ^/^^ ^, aio. 



vrt \ ''^ sin t vr cos t vr sin 3t 

— (sin t - cost) H 1 -j= _ 

4 ^ ^ 4 8 8^3 



+ 0(t 



vrt N ^ sin t vr cos t vr sin 3t 

— (sin t + cost) 1 -j= _ 

4 ^ ^ 4 8 8^3 



+ 0(t-9/8)^ 



Proof. Suppose that /(t), t > 0, is a smooth function which is bounded near t = 0. 
Let F{t) = /o /(s)ds and suppose that fait) and Fa{t) are smooth functions bounded 
near t = such that 



f{t) = fa{t)+Oit~^), F{t)= Fait) + 0(1), t^CX). 



(3.33) 



Integrating by parts the first term on the rhs of the equation 

fis)ds 



we find 



Afit) 



F(tV4) 



Afit) 



^1/4 



fis)ds 



+ 



/4 vr 



Fis)ds 



* ifis)-fais))ds , fais)ds 



Vt - tl/4 Jo 2(t-s)3/2 7,1/4 ,/t^ 

Using the estimates 



/4 ^t- s 



,1/4 



ds 



(t- 5)3/2 



or'/'), r 



(is 



/4 



0(t-^/^), t ^ oo, 



^The lower limit of integration here can be taken to be t^^^ or 0; an estimate as in (3.321 shows 
that the difference is of 0(t"^/*). 
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we conclude that 



Am 



41/4 



Fais)ds 



y^iZTtm Jo 2{t - s)3/2 7,v4 Vi 



fais)ds 



00. 



(3.34) 



Now suppose that f{t) = Hs{t) s'mt. Then equation (3.33) is satisfied with 

1 f¥ 



Fait) 



tcost 



and 



fait) = ^y|tsin(t) + ^Y^(sin(t) - sin(3t) - cos(t) + cos(3t)) 



sin 2t 3 

+ + 



3cos2t 15sin2t 



4^3/2 8^5/2 8*5/2 



8*7/2 • 



Using that 
Fit) = 



-J-tcos(t) + ^ 



1 



1 /vr 
24V 2 
cos(2t) 9sin(2t) 



4i3/2 8t3/2 



32t5/2 



9 sin(t) + sin(3t) - 3 cos(t) + cos(3t)) 



it is easy to find the asymptotics of all the terms in (3.34) up to 0(t ^/^^ except for 
the contribution from the following terms in /„: 

sin2t 3cos2t 15sin2t 

In order to handle the first of these terms, we consider an integral of the form 

sin(ns) 



/4 s'^l'^yjt — S 



ds, 



n > 0. 



We write 



sm ns 



s3/2vr 



zds 



t-1 



sm ns ) 



jl/4 s'^/'^^Jt — S 



ds + 



sin(ns) 



ds. 



(3.36) 



Integrating by parts twice, we can write the first term on the rhs as 



cos(ns) 



ns3/2^i _ g 



sin(ns) 



4s -3t 



,,1/4 n2 2s5/2(t - s)3/2 



t-1 



(1/4 



1 sin(ns) 3 (Ss^ - 12si + ht^) 

4s7/2(i_g)5/2 



ds. 



(3.37) 



Using that 8s^ — 12st + bt^ > for s G it^/^, t — 1) for t large enough and that 

3 (Ss^ - 12si + 5*2) 

7(1/4 4s7/2(t_s)5/2 



'ds = 0(t-9/S), 



we deduce that the last integral in (3.37) is of 0(t ^/^). On the other hand, the 

1 



second term on the rhs of (3.36) can be estimated as follows: 



sm ns 



it-i 4s3/2^/t^ 
In summary. 



ds 



< 



it - 1)3/2 J^_^ 4^ 



1 



.ds = on 



-3/2^ 



t 



00. 



L 



sin(ns) 
(1/4 s3/2^i — s 



ds 



cos 



(nti/4) 



n 



+ 0(t-9/8). 
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This determines the contribution to equation (3.34) from the first term in (3.35). 
An analogous argu ment shows that the two other terms in (3.35) yield contributions 
to equation (3.34) of 0{t~^^^). Equation (3.34) now gives the asymptotics of the 
function A{Hs sin) as stated in the lemma; the other functions can be analyzed in a 
similar way. □ 



Combining the known asymptotics of h, He, Hg with lemma 3.12, we find the 
following asymptotics of the last eight quantities in (3.30). 



Lemma 3.13. Let h,I,A,Hc,Hs be defined in {3.1^ and (3.20^ - (3.2^ . Then the 
following estimates are valid: 

7r3/2 (t^ 



I[hA{Hc cos)] 



t t cos 2t sin 2t 3 sin 2t 
yl"l"8 ~ 16 8 32^/3 ^ 32 
sin At cos 2t cos 2t cos 4t 1 



64^/3 32^/3 



+ 



32 



+ 



64^/3 J 



+ ki + 0{t 



^l/8x 



I[hA{Hs cos)] 



vr 



3/2 



t t cos 2t sin 2t 5 sin 2t 
V2{ 8 16 8 32^1 ^ 32 



sin4t cos2t cos2t cos4t 

H ?= + ^ + -^z:^ + 



64^/3 32^/3 32 64^/3 



+ ^2 + 0(^1/8), 



sin)] 



vr 



3/2 C ^ 



tsin2t sin2t sin2t sin4t 
V2 1 16 8 ^ 32^/3 ~ ^32^ ^ 64^/3 



+ 



cos 2t 3 cos 2t cos 4t 1 



32^/3 



32 



64\/3 J 



+ k3 + 0{t-^l 



I[hA{Hcsm)] 



7r^/2 f 3t tsin2t sin2t sin2t sin4i 

+ 



\/2 116 8 32^/3 32 64^/3 
cos 2t 5 cos 2t cos 4t 1 
^ 64 V3 J 



+ 



32\/3 



32 



+ ki + o{t-^' 



and 



I{hHc cos) 
I{hHs cos) 

I{hHs sin) 
I{hHc sin) 



+ 



vrt^ 7rt(sin 2t — cos 2t) 7rsin4t 
"16 
7rt2 



16 



64 



4Vt 



vrt 7rt(sin 2t — cos 2t) 7rcos2i 



16 16 

vr cos 4f 
^ 64 



+ 



16 



+ 



f sin t 

vrt^ 7rt(sin 2f + cos 2t) ^ 7rsin4t 



16 

-3/2>, 
+ 



16 16 ' 64 ' 4Vt 

vrt^ vrt 7rt(sin 2t + cos 2t) 7rcos2i 



16 ^ 16 

vr cos 4i 
H — + 



16 16 

2*^°^* , , /^^4-3/2^ 



64 ' 4\/t 

where {kj^\ are real constants. 



+ fes + o{r 



Proof. We consider I[hA{HcCos)\\ the other terms are treated similarly. By (3.16) 
and lemma 13.121 we have 



hA{H,cos) = f{t) + fr{t), 
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where f{t) is a known function and frit) = 0{t '^/^). Choose K > and M > such 
that 

l/rWI<^, t>K. 

Then L := fr{s)ds is a finite number and, as t — ?• oo, 

I[hA{HcCos)]- hA{HcCos)ds- f{s)ds-L = - fr{s)ds = 0{r^^^). 
Jo Jk Jt 

Thus, up to a constant, we can find the asymptotics of I[hA{HcCos)] to 0{t^^^^) by 
simply integrating f{t). This gives the result. □ 



The asymptotics of the first eight terms in (3.30) to 0{t~^^^^) now follow from 
lemma 3.11 in view of the asymptotics of lemmas |3.10| and |3.13 Indeed, all the 
functions acted on by A in these terms are at most of O(t^) as t — t- oo and lemma 



3.11 applies. 



It only remains to treat the last term in the expression ( |3.29| ) for Ihhs- This term 
is not in the form of an Abel integral. However, we can find its asymptotics using 
different arguments. 



Lemma 3.14. The last term on the rhs of (3.29) satisfies 



2 / dt'h{t' 



dt 



„cost"Hc{t") + smt"Hs{t" 



Vt-t' + 1" 

= \/27r[/i(t) cost + hit) smt + I^it) cost + h(t) s'mt] + 0{t~^^^ 
where the functions Ii, I2, 13, h are defined by 

h = dt" — / [He cos^ +{Hs + He) cos sin +Hs sm'^)du, 
Jo Vt" Jo 



(3.38) 



f cos t f 

I2 = I dt" — -1=- I {Hccos^ +{Hs — He) cossm—HsS\T?)du, 
Jo vt" Jo 

I 

{—HcCos^ +{—Hs + He) cossin+i^s sm'^)du, 



* sint" 
dt" 



h 



t" Jo 



t" Jo 



[He cos^ +{Hs + He) cos sin +Hs sm'^)du. 



(3.39a) 
(3.39b) 
(3.39c) 
(3.39d) 



Proof For the purposes of finding the asymptotics to 0{t ^/^), the function h{t') in 



implies that 
Since 



the Ihs of (3.38) may be replaced with its leading asymptotics. Indeed, lemma 3.10 
cost"He{t") + smt"Hs{t") = 0{t"), t" 00. 

1 



2 [ dt' [ dt" 
Jo Jo 



Vt-t' + 1" 



o{t-^'% 



and 



2 t dt' ^ 



t' 



dt"- 



t" 



1 t'V^ Jo Vt-t' + 1 



0{t 



-l/2^ 



it follows that to 0{t ^^'^) the asymptotics of the Ihs of (3.38) coincide with those of 
the following expression: 



2 / dt\/-(sin(t') + cos(t')) / dt 



„ cos t" He{t") + smt" Hsit"^ 



Vt-t' + 1" 
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2 dt' ^ 



*' ^^„mst"Hc{t") + smt"Hs{t") 



1 2t'3/2 7q - Vt-t' + 1" 

Changing the order of integration and performing the integral with respect to dt\ we 
can write this expression as 

27r / dt"[cost"Hc{t")+s\nt"Hs{t")\ 
Jo 

X I (sin(t + t") + cos(t + t")) (^C ^) - ^ (/f ^ 
+ (sin(t + f) - cos(t + t")) (^s(^^Vt^ - s{J^Vt^ 

+2 f dt" [cos t"Hc{t") + sin t"Hsit")] 
Jo 

-2 [ dt"[cost"Hc{t")+smt"H,{t")]j^^—j;^^y=. (3.40) 



t + t" 
t-t" 



{t + t")Vtt" 



The first integral in (3.40) becomes after integration by parts 

V^TT J dt" du[cos uHc{u) + sin uHs{u)]{sm{t + u) + cos{t + u)) 



cos t" 



\Jo 



+ V2tt / dt" I / du[cos uHc{u) + sin uHs{u)]{sm{t + u) — cos{t + u)) 



t" 

sint" 



t" 



After expanding the above trigonometric functions it follows that this expression 
equals 

V27r [h (t) cos t + I2 (t) sin t + h (t) cos t + h{t) sin t] . 

It therefore only remains to prove that the last two integrals in (3.40) are of 0(t~^/^). 
For the next to last integral, this is clear. To see that the final integral is of 0(t^^/^), 
we use repeated integration by parts to write it as follows: 

/■*" t-t" * 

— 2 / [cosuHc{u) + sinuHs{u)]du 



{t + t")Vtt" 



t"=i 



t" rr 



+ 2 / / [cos uHc{u) + sin uHs{u)]dudT-— 



d f t-t" 



1 Ji 

ft ft" fT 



dt" \(t + t")Vti" 



[cos uHc{u) + sinuHs{u)]dudT - 



t - t 



t"=i 



dt"'^ \ (t^t")^fW' 



(3.41) 



The following facts imply that the last line of (3.41) is of 0{t ^/^): 
• /i* ^^[cos ullciu) + svnuH s(u)\dudT is of 0(t") as i" — )• 00, 



'^^ ^ ^ > for aU < t" < t, 

0(t-i/2). 



''i J+//J.// li?' I t~t 



Jl dt"t" 



dt" 



{t+t")Vtt" 



Consequently, the expression in (3.41) equals 

ft rr 



1 



1 Jl 



[cosuHciu) + sinuHs{u)]dudT + 0{t ^/^) = 0(t 



-l/2^ 



□ 
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We can no w com plete the proof of the theorem. Using the expressio ns for He and 
the asymptotics of the integrals {Ij}i defined in (3.39) can be 



3.10 



Hg of lemma 

computed up to a constant. This yields the following formulas: 



h 
h 
h 
h 



TT 

32 \ 

TT 

32 



f x/tsint \/tsm3t 3\/tcost 
<r' sint-—: — h 



-Vtsint 



2 6 2 

^/i cos t \/t cos 3t 



+ A;9 + 0(t-i/2), 



6 



vr f v/isint \/tsin3t /- i , ^/ i/9n 

2 D 



32 
vr 
32 



(3.42) 



t^/^ cos t + 



3\/tsint \/tcost \/tcos3t 



+ 



6 



+ A:i2 + 0(f 



-1/2N 



where {fcjjg^ are real constants. 

Combining all of the above, we arrive at the following asymptotics of Tj, j 
1, . . . , 7, as t — )• oo: 



Ti = |9(cost - sint) - \/3(cos3t - sin3t)| + 0(t"^ 



Ta = ^'-^|6sint + (\/3-3)sin3t-12cost- (\/3 - 3) cos 3t| +0{r^/^), 
48 I- J 



48 
sint 

1 -i 



(2t - sin 2t - cos 2t) + 0(t"^/^), 
|4i7rt(sint — cost) 
+ 2iiT ^3 sin t + cos t + 

(\/3- l)(cos(3t) -sin3t))} + 0(t-5/2)^ 

Ts = ^7^|-47rt(cost - 3sint) + 32^{{ki + /ca) cost + (A;2 - fc4)sint) 

+ sint(7r — 32(A;5 + k^ + k-j — kg)) + cost(7r — 32(A;5 — /cg + A;? + kg)) 
+ 2 (^V3 - 7r(cos3t - sin3t)| + 0{r^/^^), 

AiTTt{3 sin t + cos t) + 4i sin t (S/cs + 8A;6 - Skj + 8/^8 - 27r) 



Tf 



647r 



1 -i 
647r 



1 -z 



647r 



+ 4i cos t(8A:5 - 8{k(i + kj + kg)) + 2i (^\/3 - 7r(cos 3t - sin 3t) | + ©(t^^/^*^), 

Ty = — — <^ 47rt(sint + cost) + sint(32A;io + 32A;i2 - GAk^ - GAky - vr) 
647r L 

+ cost(32A;ii - GAk^ - Ukj + 32/09 - vr) 

+ 2 {y3 - 1) 7r(cos3t - sin3t)| + 0{t-^/^^). 

Substituting these expressions into ( 3.17| ), we find that 5(13 satisfies (3.14c) with ci 
and C2 given by 



ci 



C2 



l-i 
2tt 



l-i 
2tt 



{ki + fcs) - (3 - i)k^ + (1 - i)kQ - (3 + i)A:7 



(1 + i)kg + fcg + fell + ( ^ + ^ I TT 



(^2 — ki) — (3 — i)k^ — (1 — — (3 + i)kj 



(3.43a) 
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+ {l + i)k8 + kio + k 



Cl2 



1 lli 

- H I vr 

8 16 



(3.43b) 



where {/cj}} are defined in lemma 3.13 and equation (3.42 ). This completes the proof 



of theorem 13.2 



□ 



Remark 3.15. 1. Neither of the terms T3-T7 is asymptotically periodic by itself, 
however their sum is. 

2. The leading asymptotic behavior of 0{t) of T^, Tq, and T7 can be obtained 
by simply replacing h{t) by its leading asymptotics ^y^{cost + sint) in (3.18) and 
computing all the integrals in terms of C{z) and S{z). 

Appendix A. An alternative system 



Eliminating the unknown functions gi and gi from the system (2.3), we find after 
simplification the following equations: 

/3M| + ( L2 + ^5oAfi ) [Mu - Mis - 250^2] + ^^50^2 [Mat + M2, - 2A5o^i] 



M2[Lu-Lu] 



and 



A/3M2 + f Li - ) [M2t + M2s - 2Ago^i] 



goMi[Mit - Mu - 2goL2] 



2 -J' 2 

= Mi[L2t + L2s\. 

Replacing in the above Li and L2 using (2.14a) and (2.14b), we find 

Mi{t,T)dT 

-t y/s-T 



M2{dt - ds)ds 



{Mu - Mu)ds 

Mi{t,T)dT 



M2{t,T)dT In 



2 



M2it,T)dT 
-t \/s-T 







and 



' M2{t,s)dT 
t y/s-T 

a 1 50 1' fM2a 



(M2t + M2.)a 

^ Mi{t,s)dT 



s — r 



Mia, 



t Vs-T ' "V 2 
' M2{t,s)dT 



UMt 



-t \/s-T 
2 



Mi{t,s)dT 

-t VS-T 



0, 



This provides an alternative nonlinear system for Mi and M2. However, it is not 
clear if this system is perturbatively solvable. 



Appendix B. Derivation of (2.7) 



In order to solve equation (|2.6|) we introduce the characteristic coordinates ^, and 



t+s 
2 ' 



V 



S = ^ - 7?. 
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dt-ds = dr,, 

The restriction —t<s<t imphes that 

0<t- s <2t, 

thus 
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Equation (2.6) yields 
Hence, 



< r? < T 

dF 
drj 



dt + ds = d^. 
0<t + s <2t, 
< ^ < T. 



F{t, s) = F{C + Vp,C- Vp) + r /(e + e - v)dv', 0<Vp<T. 

J rip 



Thus, 

F{t,s) 



F 



t + s 



t + s 



rip + 



t-3 

2 



f 



Vp 



t 



2 2 

Letting r/p = and employing the change of variables 

t + s 



— + r]', ^ \ ^ - rj' ] drj'. (B.l) 



+ rj 



equation (B.l) becomes equation (2.7). 



Alternatively, letting rjp = and employing the change of variables 

t + s 



V 



equation ( |B.l ) yields 

't + s t + s 



Similar considerations are valid for equation (2.9). 



-t <s <t, <t<T. 



Note added: After this paper was written we were pleased to find out about [T], 
where the equation obtained from the KdV by adding the term aq{x,t), a > 0, 
representing damping is considered. For this equation, it is shown via techniques 
of functional analysis, that if the small amplitude boundary forcing is periodic of 
period T, then the solution q{x, t) eventually also becomes periodic of period T at 
each spatial point. 
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